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^ Abstract 

O ' A new formalism of beam-optics and polarization has been recently pre- 

c/2 . sented, based on an exact matrix representation of the Maxwell equations. 

! This is described in Part-I and Part-II. In this Part, we present the applica- 

j>^| tion of the above formalism to the specific example of the axially symmetric 

I graded index fiber. This formalism leads to the wavelength-dependent mod- 
ifications of the six aberrations present in the traditional prescriptions and 
further gives rise to the remaining three aberrations permitted by the axial 

> . symmetry. Besides, it also gives rise to a wavelength-dependent image rota- 

' tion. The three extra aberrations and the image rotation are not found in 

00 

I any of the traditional approaches. 

vn ' 
o 

'. 1 Introduction 

o 

c/3 ' In Part-I and Part-II we presented the exact matrix representation of the 

Maxwell equations in a medium with varying permittivity and permeabil- 
CL|! ity m, 0]. From this we derived an exact optical Hamiltonian, which was 

shown to be in close algebraic analogy with the Dirac equation. This enabled 
us to apply the machinery of the Foldy-Wouthuysen transformation and we 



' obtained an expansion for the beam-optical Hamiltonian which works to all 



orders. Formal expressions were obtained for the paraxial Hamiltonian and 
the leading order aberrating Hamiltonian, without assuming any form for 
the refractive index. Even at the paraxial level the wavelength-dependent 
effects manifest by the presence of a matrix term coupled to the logarithmic 
gradient of the refractive index. This matrix term is very similar to the spin 



term in the Dirac equation and we call it as the polarizing term in our formal- 
ism. The aberrating Hamiltonian contains numerous wavelength-dependent 
terms in two guises: One of these is the explicit wavelength-dependent terms 
coming from the commutators inbuilt in the formalism with A playing the 
role played by h in quantum mechanics. The other set arises from the the 
polarizing term. 

Now, we apply the formalism to specific examples. One is the medium 
with constant refractive index. This is perhaps the only problem which can 
be solved exactly in a closed form expression. This is just to illustrate how 
the aberration expansion in our formalism can be summed to give the familiar 
exact result. 

The next example is that of the axially symmetric graded index medium. 
This example enables us to demonstrate the power of the formalism, repro- 
ducing the familiar results from the traditional approaches and further giving 
rise to new results, dependent on the wavelength. 



2 Medium with Constant Refractive Index 



Constant refractive index is the simplest possible system. In our formalism, 
this is perhaps the only case where it is possible to do an exact diagonal- 
ization. This is very similar to the exact diagonahzation of the free Dirac 
Hamiltonian. From the experience of the Dirac theory we know that there 
are hardly any situations where one can do the exact diagonahzation. One 
necessarily has to resort to some approximate diagonahzation procedure. 
The Foldy-Wouthuysen transformation scheme provides the most convenient 
and accurate diagonahzation to any desired degree of accuracy. So we have 
adopted the Foldy-Wouthuysen scheme in our formalism. 

For a medium with constant refractive index, n (r) = ric, we have. 



He = -ricP + i {MyPj: - M^py) , 



(1) 



which is exactly diagonalized by the following transform. 



T 



exp [i {±if3) Oe 

exp [^i/3 {Myp^ - M,py) 9] 




sinh 9) 



(2) 
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We choose, 



then 



tanh(2|p^|^) = ^ (3) 

Tin 



l2P, {n, + P,) 



where — +^ (nl — p^) . Then we obtain, 

pjdiagonal ,T~ 

= T+{-n,P + i{Myp,-M,py)}T- 

= -{nl-pl}"^P (5) 

We next, compare the exact result thus obtained with the approximate ones, 
obtained through the systematic series procedure we have developed. 
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j^diagonal 

Knowing the Hamiltonian, we can compute the transfer maps. The trans- 
fer operator between any pair of points {{z",z') \z" > z'} on the z-axis, is 
formally given by 

\i^{z",z')\=fiz",z')mz",z')) , (7) 

with 

iA^r(z", z') = nf{z", z') , f{z\ z') = x, 
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T(z", z') = pi exp 



dzH{z) 



+ ... , 



dzTii^z) 



dz / dz'n{z)n{z') 



(8) 



where X is the identity operator and p denotes the path-ordered exponential. 
There is no closed form expression for 'T{z" ^ z') for an arbitrary choice of the 
refractive index n{r). In such a situation the most convenient form of the 
expression for the 2;-evolution operator T{z",z'), or the z-propagator, is 



f{z",z' 



exp 



(9) 



with 



f{z", z' 



dz7i{z) 
+ ... , 



dz / dz' n{z),n{z') 



(10) 



as given by the Magnus formula §]. We shall be needing these expressions 
in the next example where the refractive index is not a constant. 

Using the procedure outlined above we compute the transfer operator, 



Uc (^^out) ^in) — GXp 



where, Az 



exp 



[2^0Ut) Ziyi 



+jn,Az { 1 



X 



2nl 8 



nt 



(11) 



I. Using ([TT|), we compute the transfer maps 



out 



1 





1 



■Az 



(ri.) 
iPi.) 



(12) 



The beam-optical Hamiltonian is intrinsically aberrating. Even for simplest 
situation of a constant refractive index, we have aberrations to all orders. 
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3 Axially Symmetric Graded Index Medium 

The refractive index of an axially symmetric graded-index material can be 
most generally described by the following polynomial (see, pp. 117 in 0]) 

n (r) = no + a2{z)r\ + a4{z)r'[ -\ , (13) 

where, we have assumed the axis of symmetry to coincide with the optic-axis, 
namely the 2;-axis without any loss of generality. We note, 

£ = - {a2iz)rl + ai{z)r\_ + ■■■,} P - iA/3S ■ u 

6 = i{Myp^-M^py) 

= P{M^-p^) (14) 

where 

S-n = -— a2(^)S^-rx--^(-^a2(^))S,ri (15) 
riQ 2nQ \dz J 

To simplify the formal expression for the beam-optical Hamiltonian Ti.^^^ 
given in (24-25) in Part-II, we make use of the following: 

(M^-pJ = pI, d' = -pl, §;O = 0, 
(M^-K)^i(Mx-Px) = i(ripl + plri)+2A/3L, + 2A2, (16) 

where, is the angular momentum. Finally, the beam-optical Hamiltonian 
to order (^-^p]_^ is 

Ho,p = -no + ^p']_ - a2{z)r']_ 
Zno 

1 -4 



^0.(4) = ^Pl. 



[rlPl+plr] 
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(A) 
0,(2) 



" 



A 



2nl 



^a2{z)L, + —al{z)rl 



H. 



(A) 
0,(4) 



3 1 J.""2(^) } P^z 



iAM , f 



2r;,Q [ dz 
iA 



a2iz)P S_L ■ r±,p]_ 



iA f d 



[ dz 



a2iz) \ (3J:, 



2 - 2 



+ 



(17) 



where -B] = {AB + SA) and '■ ■ ■' are the numerous other terms arising 
from the polarization term. We have retained only the leading order of 
such terms above for an illustration. All these matrix terms, related to the 
polarization will be addressed elsewhere. 

The reasons for partitioning the beam-optical Hamiltonian 7i in the above 
manner are as follows. The paraxial Hamiltonian, Ho,p, describes the ideal 
behaviour. Ho, (4) is responsible for the third-order aberrations. Both of these 
Hamiltonians are modified by the wavelength-dependent contributions given 



in Hq and respectively. Lastly, we have H^^''^\ which is associated 

with the polarization. 

From these sub-Hamiltonians we make several observations: The term 
-a2{z)Lz which contributes to the paraxial Hamiltonian, gives rise to an 



2n: 



image rotation by an angle 0{z): 



dza2{z) 



(18) 



This image rotation (which need not be small) has no analogue in the square- 
root approach |^ and the scalar approach |^ . 
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The Hamiltonian Ho,{a) is the one we have in the traditional prescriptions 

and is responsible for the six aberrations. -^o'*'(4) J^odifies the above six aber- 
rations by wavelength-dependent contributions and further gives rise to the 
remaining three aberrations permitted by the axial symmetry. Before pro- 
ceeding further we enumerate all the nine aberrations permitted by the axial 
symmetry. The axial symmetry permits exactly nine third-order aberrations 
which are: 



Symbol 

C 

K 

k 

A 

a 

F 

D 

d 

E 



Polynomial 
{p± -r^ + rj 



(Plr 

^4 



r_L + r 

2 I „2 



■ P±) Lz 



Name 

Spherical Aberration 
Coma 

Anisotropic Coma 
Astigmatism 
Anisotropic Astigmatism 
Curvature of Field 
Distortion 

Anisotropic Distortion 
Nameless? or POCUS 



The name POCUS is used in on page 137. 

The axial symmetry allows only the terms (in the Hamiltonian) which 
are produced out of, p^, r^, ■ + rj_ -p^) and L^. Combinatorially, 
to fourth-order one would get ten terms including L^. We have listed nine of 
them in the table above. The tenth one namely, 



Li 



1 



2 2 



2 -2 



4 



(19) 



So, Ll is not listed separately. Hence, we have only nine third-order aberra- 
tions permitted by axial symmetry, as stated earlier. 
The paraxial transfer maps are given by 




out 




(20) 



where P, Q, R and S are the solutions of the paraxial Hamiltonian (|17|) . The 
symplecticity condition tells us that PS — QR =1. In this particular case 
from the structure of the paraxial equations we can further conclude that: 
R = P' and S = Q' where ' denotes the z-derivative. 
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The transfer operator is most accurately expressed in terms of the paraxial 
solutions, P, Q, R and S, via the interaction picture M. 



f {z , Zq) 



exp 
exp 



A 



f {z , Zo) 



'-l[c{z",z')pi 

+K {z" , z') [p\ , (K -r^ + r^- p^) 

+k{z" ,z')PlL, 
+A {z" , /) ■ + • p^f 
+a {z" , z ) (p_L ■ r_L + r_L ■ ^2 
+F{z\z')[plrl + rlpl) 

+D {z" , 2') [rl , (p^ ■ + ■ K) 

+d {z" , /) rlL, 

+E iz" , z') rl 



The nine aberration coefficients are given by, 

1 „/! a2{z 



C iz" , z') 



dz 



,S'-^Q'S'~a,{z)Q' 







a2{z)a4{z)Q^ 



K iz" , z' 



(21) 



2ng 

^RS'^ - '^^QS{PS + gi?) - ai{z)PQ' 



2nl 



a2{z)ai{z)PQ'^ 



k {z" , z') 
A {z" , 



A 



dzal{z)Q'^ 



dz 



2 c2 



R'S 



2n2 



PQRS - a4{z)P^Q'^ 



a2{z)a4{z)P^Q^ 



2nl 
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a {z" , z') 
F {z" , z') 



D [z" , z') 



d {z" , z') 
E {z" , z') 



Thus we see that the current approach gives rise to all the nine permissible 
aberrations. The six aberrations, familiar from the traditional prescriptions 
get modified by the wavelength-dependent contributions. The extra three 
(fc, a and d are all anisotropic!) are all pure wavelength-dependent aber- 
rations and totally absent in the traditional square-root approach 0, Q and 
the recently developed scalar approach P, |^. A detailed account on the 
classification of aberrations is available in |[9|]-[p!2|. 

4 Conclusions 

In Part-I and Part-II, we developed an exact matrix representation of the 
Maxwell equations which became the basis for an exact formalism of Maxwell 
optics. An exact optical Hamiltonian, with an algebraic structure in direct 
correspondence with the Dirac equation of the electron was derived. Then 
following a Foldy-Wouthuysen transformation technique, a procedure was 
developed to obtain the beam optical Hamiltonians to any desired degree 
of accuracy. Formal expressions were obtained for the paraxial and leading 
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A 

2^ 



A 

2^ 



dzal{z)PQ 



dzl—^R'S^ - ^^{P^S' + Q'R^) - a,iz)P^Q'' 

A^ 1 

+Y^o:2{z)a4{z)P^Q^ I 



dz { -^R^S - -i\^PR(PS + QR) - aAz)P^Q 

A^ 1 

Ta2{z)ai{z)P^Q \ 



2ni 



dzal{z)P^ 



dz I -1-R^ - ^^P^R'' - a,{z)P^ 
\8nl 2nl ^ ' 

+^«2(^)a4(^)p4 . (22) 



order aberrating Hamiltonians, without making any assumption on the form 
of the refractive index. In this Part we look at the apphcations of the above 
formahsm. 

First of the two examples is the medium with a constant refractive index. 
This is perhaps the only problem which can be solved exactly, in a closed 
form expression. This example is primarily for illustrating certain aspects of 
the machinery we have used. 

The second, and the more interesting example is that of the axially sym- 
metric graded index medium. For this example, in the traditional approaches 
one gets only six aberrations. In our formalism we get all the nine aberra- 
tions permitted by the axial symmetry. The six aberration coefficients of the 
traditional approaches get modified by the wavelength-dependent contribu- 
tions. 

It is very interesting to note that apart from the wavelength-dependent 
modifications of the aberrations, this approach also gives rise to the image 
rotation. This image rotation is proportional to the wavelength and we have 
derived an explicit relationship for the angle in (|18[). Such, an image rotation 
has no analogue/counterpart in any of the traditional prescriptions. It would 
be worthwhile to experimentally look for the predicted image rotation. The 
existence of the nine aberrations and image rotation are well-known in axi- 
ally symmetric magnetic electron lenses, even when treated classically. The 
quantum treatment of the same system leads to the wavelength-dependent 
modifications [l^ . 

The optical Hamiltonian has two components: Beam-Optics and Polar- 
ization. We have addressed the former in some detail and shall do the later 
soon. The formalism initiated in this article provides a natural framework 
for the study of light polarization. This would provide a unified treatment 
for the beam-optics and the polarization. It also promises a possible gen- 
eralization of the substitution result in ||T^. We shall present this approach 
soon 

The close analogy between geometrical optics and charged-particle has 
been known for too long a time. Until recently it was possible to see this 
analogy only between the geometrical optics and classical prescriptions of 
charge-particle optics. A quantum theory of charged-particle optics was pre- 



sented in recent years JTS], ^ |I3]. With the current development of the 
non-traditional prescriptions of Helmholtz optics P, |^ and the matrix for- 
mulation of Maxwell optics (in these three Parts), using the rich algebraic 
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machinery of quantum mechanics it is now possible to see a parallel of the 
analogy at each level. The non-traditional prescription of the Helmholtz op- 
tics is in close analogy with the quantum theory of charged-particles based on 
the Klein-Gordon equation. The matrix formulation of Maxwell optics pre- 
sented here is in close analogy with the quantum theory of charged-particles 
based on the Dirac equation. We shall narrate and examine the parallel of 
these analogies soon [pT|. 
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